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Abstract. 

We present SWS and LWS ISO observations towards a strip 
across the photodissociation region (PDR) in the reflection 
nebula NGC 7023. SWS02 and LWSOl spectra have been taken 
towards the star and the brightest infrared filaments located 
NW and SW from the star (hereafter NW and SW PDRs). In 
addition, SWS02 spectra have been taken towards two interme- 
diate positions (NWl and SWl). This has provided important 
information about the extent and spatial distribution of the 
warm Hg and of the atomic species in this prototypical reflec- 
tion nebula. 

Strong emission of the [Sill] 34.8 pm line is detected to- 
wards the star. While aU the PDR tracers (the [CII] 157.7 /xm, 
[OI] 63.2 and 145.6 /xm, [HI] 21cm and the H2 rotational lines) 
present a ring-like morphology with the peaks toward the NW 
and SW PDRs and a rninirrmrn around the star, the Sill omis- 
sion is filling the hole of this ring with the peak towards the 
star. This morphology can only be explained if the Sill omis- 
sion arise in the lowest extinction layers of the PDR (A„ < 2 
mag) and the HII region. At least 20% - 30% of the Si must be 
in gas phase in these layers. For A„ > 2 mag, the Si is mainly 
in solid form {5 Si — -1.3). 

The NW and SW PDRs have very similar excitation con- 
ditions, high density filaments (n~ 10® cm~^) immersed in a 
more diffuse interfllament medium (n~ 10* cm~^). In both, the 
NW and SW PDRs, the intensities of the H2 rotational lines 
can only be fitted by assuming an ortho-to-para-H2 ratio lower 
than 3 in gas with rotation temperatures from 400 to 700 K. 
Therefore, there is a non-equilibrium OTP ratio in the region. 
Furthermore, the comparison between the OTP ratio derived 
from H2 vibrational lines and the pure H2 rotational lines, 
shows that the OTP ratio increases from ~ 1.5 to 3 across the 
photodissociation region with larger values in the loss shielded 
gas (A^< 0.7 mag). This behavior is interpreted as a conse- 
quence of an advancing photodissociation front. 

We have not detected the OH, CH and CH2 lines towards 
the observed positions. This is consistent with the weakness of 
these lines in other sources and can be explained as a conse- 
quence of the small beam filling factor of the dense gas in the 
LWS aperture. The CO J=17— »16 line has been tentatively 
detected towards the star. 
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1. Introduction 

NGC 7023 is a prototypical photodissociation region that has 
been largely studied in the last 10 years. It is a reflection nebula 
niuminated by the Herbig B3Ve star HD 200775. HD 200775 is 
located in a cavity of the molecular cloud whose sharp edges de- 
lineate perfectly the optical nebula (Fuente et al. 1992, Rogers 
et al. 1995, Fuente ot al. 1998). Observations of CI and Oil 
by Chokshi et al. (1988) reveal that a dense PDR is formed 
in the walls of this cavity with its peak located 50" NW from 
the star. Chokshi et al. (1988) carried out the first model for 
this region and estimated a UV field of Go(A> 912 A)=2.4 10^ 
in units of the Habing field and a density of n = 10'' cm"''. 
Because of its edge-on geometry and proximity (d ~ 440 pc), 
this PDR turned out to be one of the best sites to study the 
physical and chemical processes taking place in a PDR. Fig. 1 
shows the column density map of HI as derived from the [HI] 
21 cm line superposed on the integrated intensity map of the 
^^CO J=1^0 line. The HI peak is shifted « 20 " relative to 
the molecular ridge, showing the layered structure expected in 
a PDR (Fuente et al. 1998). 

Very high angular resolution images of the region in the 
Extended Red Emission (ERE) and vibrationally excited H2 
lines show that the PDR has a filamentary structure with very 
bright thin filaments located 50" NW and 70" SW from the star 
(Sellgren et al. 1992, Lemaire et al. 1996). Hereafter, we will 
refer to these positions as NW PDR and SW PDR respectively. 
The NW PDR is spatially coincident with the peak of the CI 
and on lines and has been extensively studied in atomic and 
molecular linos (Chokshi et al. 1988; Fuente et al. 1993, 1996, 
1997; Lemaire et al. 1996,1999; Martini et al. 1997; Gerin et al. 
1998). Fuente et al. (1993, 1997) observed that the CN/HCN 
and CO"'"/HCO"'" abundance ratios increases by a factor 10 and 
> 100 respectively towards this position relative to the values 
in the molecular cloud. Both, the CN/HCN and CO+/HCO+ 
abundance ratios are expected to increase in the lower extinc- 
tion layers of PDRs (Fuente et al. 1993, Sternberg & Dalgarno 
1995). An interferometric image of the HCO"'' line towards the 



NW PDR show that the PDR has also a filamentary structure 
in molecular emission with high density filaments (n > 10' 
cm^"*) embedded in a more diffuse medium (n ~ 10^ cm^'^) 
(Fuente et al., 1996). The HCO+ filaments are spatially coin- 
cident with the filaments seen in the neax-infi'ared continuum 
images and in the emission of the vibrational H2 lines. 

The SW PDR is fainter by a factor 2-3 than the northern 
one and has been less studied in both, molecular and atomic 
lines. Gerin et al. (1998) reported observations of this PDR in 
CO and CI fines. 

We have observed a strip which joins the NW PDR, the star 
and SW PDR using the SWS and LWS instrument on board 
of ISO. These observations have provided a very important 
information about the extent and spatial distribution of the 
warm H2 and the atomic species within the PDR. 

2. Observations 

The observations were made using the Long- Wavelength Spec- 
trometer (LWS) (Clcgg ct al. 1996, Swinyard et al. 1996) and 
the Short- Wavelength Spectrometer (SWS) (de Graauw et al. 
1996) on board the Infrared Space Observatory (ISO) (Kessler 
et al. 1996). The LWS observations were carried out in full grat- 
ing scan mode (LWSOl AOT) which provides coverage of the 
43 - 90.5 fiin range with a spectral resolution of 0.29 fim and 
of the 90.5 - 196.5 /xm range with a resolution of 0.60 /im with 
ten independent detectors. Since the instrumental beam size is 
74" - 90", we have observed only three positions of the nebula, 
the NW PDR, the star and the SW PDR. In Table 1, we give 
the coordinates of the observed positions. The data were pro- 
cessed using version 6 of the off-line pipeline (OLP V6). The 
uncertainty in the calibration is of about 30 % (Swinyard et 
al. 1996) 

We have observed one HI recombination line (Bra), five H2 
pure rotational lines (S(0),S(1),S(2),S(3), S(4) and S(5) v=0 0) 
and four ionic and atomic fine structure lines ([SIV] 10.5 /xm, 
[Fell] 26.0 fxia, [SIII] 33.5 fim and [Sill] 34.8 /xm) using SWS 
AOT02. This mode provides a spectral resolution of A/ A A ~ 
1000-2000. All the observed lines are unresolved at this spectral 
resolution. The aperture of the instrument in those line ranges 
from 14" X 20" to 20" x 33". In order to study the spatial 
distribution of these lines we have observed 5 positions along 
the strip that joins the NW PDR, the star and the SW PDR. 
These positions are marked in Fig. 1, the observed spectra are 
shown in Fig. 2 and absolute coordinates arc given in Table 1. 
Data reduction are carried out with version 7.0 of the Off Line 
Processing routines and the SWS Interactive Analysis at the 
ISO Spectrometer Data Center at the Max-Planck Institut fiir 
Extraterrestrische Physik. The uncertainties in the calibration 
are of about 20 % at short wavelengths (2 - 20 /im) and of 
about 30% at longer wavelengths (Salama et al. 1997). 

3. Results 

The LWS and SWS observational results are shown in Table 2 
and 3. The spectra of the SWS observations are shown in Fig. 
2. In order to study the spatial distribution of the different 
linos, wo have plotted in Fig. 3 the normalized intensities of 
the [HI] 21cm line (Fuente et al. 1996,1998), the H2 S(l) v=0- 
and [Sill] 34.8 /xm lines as a function of the distance from 
the star. 



3.1. Recombination lines 

The Bra line has only been detected towards the star. Con- 
tinuum emission at 3.6 cm and 6 cm was detected by Skinner 
et al. (1993). They interpreted the emission as arising in an 
ionized anisotropic wind. Later, Nisini et al (1995) observed 
the Pa/3, Br7, Bra, Pf7 and Pf/3 recombination lines with an 
aperture of 10" and modeled the wind with a mass loss rate of 
3.9 10"'^ Mo yr"^ and a terminal velocity of 280 kms"^ We 
have measured an intensity of the Bra line of 1.4 10~^ erg s~^ 
cm~^ sr^^ which agrees within the calibration uncertainties 
with the data of Nisini et al. (1995). 

3.2. The CII and 01 lines 

The [CII] 157.7 ^im and the [01] 63 and 145.6 ^m lines have 
been detected towards all the observed positions. The intensi- 
ties and line intensity ratios are shown in Table 2. The emis- 
sion of the three lines is maximum in the NW PDR and de- 
creases towards the south, being the SW PDR fainter by a 
factor of 2 than the NW PDR. Along the observed strip, the 
OI(63^im)/OI(145.6^im) and OI(63^m)/CII(157.7 ^m) ratios 
are uniform and take the values ~ 10 and ^ 2 respectively. 
This means that at the scale of the LWS observations, the 
physical conditions of the nebula are quite uniform. High an- 
gular resolution observations show that this PDR is formed 
by high density filaments immersed in a more diffuse medium 
(Sellgren et al. 1992, Lomaire et al. 1996, Fuente et al. 1996). 
Because of the low angular resolution of the LWS observations, 
they are very likely tracing the more extended diffuse medium. 

3.3. Other ionic and atomic fine structure lines: Sill 

The spatial distribution of the Sill line is very different from 
that of the CII, OI and HI lines (see Fig. 3). While the CII, OI 
and HI lines peak towards the NW PDR, the Sill line peaks 
towards the star, with the intensity towards the star ^ 3 times 
larger than towards the NW PDR. This different spatial dis- 
tribution cannot be explained by the different angular resolu- 
tion of the CII and 01 observations. The [HI] 21cm line map 
which has an angular resolution similar to that of the Sill line 
presents a ring-like morphology with the peaks towards the 
NW and SW PDRs and a minimum towards the star, just in 
the peak of the Sill line. The implications of the peculiar spa- 
tial distribution of the Sill emission are discussed in Section 
4.1. We have not detected the [SIII] 33.480 /xm, [SIV] 10.51 
Hm and [Fell] 25.9 fj,m lines at any position (see Table 3 for 
upper limits). 

3.4. H2 rotational lines 

We have observed the v=0-0 H2 rotational lines from S(0) to 
S(5) towards the positions shown in Fig. 1. The intensities of 
the H2 lines for all positions arc shown in Table 2. Towards 
the positions NWl, SWl and the star, wo have detected only 
the S(l) line which prevents from any excitation analysis. The 
most intense H2 emission is found towards the NW and SW 
PDRs. These positions are spatially coincident with the fila- 
ments observed in H2 fluorescent emission which are located 
in the interface between the atomic and the molecular gas. 
The H2 emission is much weaker towards the HI clump (posi- 
tion NWl) and the SWl position, suggesting that this region 
is mainly atomic. In contrast with the spatial distribution of 
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Table 1. Observations 



Position Coordinates (2000) 



NW PDR 


21''01™32=.5 


68°10'27".5 


NWl 


2l''01'°34".5 


68°10'05".0 


Star 


2l'^01'"36=.9 


68°09'48".3 


SWl 


2l''01'"36=.9 


68° 09' 10". 


SW PDR 


2l''01'°32".3 


68°08'45".0 



iil which is more intense towards NWl than towards SWl, 
the spatial distribution of the H2 emission is quite symmetric 
around the star. 

3.5. Molecular lines 

Several molecular lines of species that are expected to be abun- 
dant in PDRs (CO,OH,CH,CH2,H20) have observable lines in 
the wavelength range covered by the LWSOl spectrum. In the 
LWSOl spectrum towards NW PDR we have tentatively de- 
tected a weak line centered at 153.927 jira with a S/N ratio of 
~ 3. A weak line centered at 153.637 /xm with a S/N ratio of 
5 appears also in the spectrum towards the stax. Both spectra 
are shown in Fig. 4. Taking into account that the spectral res- 
olution of the LWS at these frequencies is 0.6 /xm, both lines 
can be identified with the CO ,J = 17^16 line (153.267 /im). 
However, since other CO lines have not been detected towards 
these positions, we cannot discard the possibility of a misiden- 
tification or a spurious detection. Towards the SW PDR, we 
have not detected any molecular line with typical upper limits 
of ~ 5 10~^^ - 5 10"^'' erg s~^ cm"'^ ^m^^ depending on the 
detector being the best limits for detectors 8 to 9 (161.2 to 
196.7 urn). 

4. Discussion 

The NW PDR was modeled for the first time by Chokslii ct al. 
(1988) with an incident UV field of Go=2.4 10^ and a density of 
n = 10* cm~^. Models of this region based on H2 ro- vibrational 
lines (Martini et al. 1997), CO, CI and CII lines (Gerin et al. 
1998) and H2 rotational lines (Fuente et al. 1999) give slightly 
higher values of the UV field, Go= 10'*, and density, n — 10^ 
- 10^ cm^'\ We can estimate the UV field and density of the 
PDR based on the observed 01 and CII lines. Along the ob- 
served strip, the 01(63 /Lim)/OI(145.6 /xm) ratio is 10. This 
value is lower by a factor ~ 2-3 than those predicted by PDR 
models for a region with a UV field of 10'^ and a density of 
lO'* cm"^ (see e.g. Kaufman et al. 1999, Burton et al. 1992 
and references therein). Low 01(63 /jm)/OI(145.6/um) ratios 
have been observed towards other PDRs and interpreted as 
the consequence of the self-absorption of the 63 /im line by 
the foreground gas (Lorenzetti et al. 1999). This interpreta- 
tion is not plausible for our region in which the extinction by 
the foreground dust is < 2 mag (Fuente et al. 1998, Gerin et 
al. 1998). We have no explanation for the low observed values 
of the 01(63 Mm)/OI(145.6Atm) ratio. 

Because the [OI] 63 /xm and [CII] 157.7 /im lines have differ- 
ent critical densities, the OI(63/im)/CII(157.7^m) ratio con- 
stitute a good density tracer for densities 10^ - l(f cm~^. The 
value of the 01(63 /im)/CII(157.7/im) is ~ 2 towards all the 
observed positions. This value is consistent with model pre- 



dictions for Go ~ 10^ and n ~ 10^ - 10* cm which are the 
physical parameters derived by Chokshi et al. (1988) from the 
same lines and better angular resolution observations. Den- 
sities around 10* cm^'^ are also derived from the rotational 
lines of CO and its isotopes (Fuente et al. 1993, Rogers et al. 
1995, Gerin et al. 1998) and from the VLA HI observations 
(Fuente et al. 1993,1998). However, densities as high as n = 
10® cm~^ are derived from the H2 ro-vibrational and rota/- 
tional lines (Lemaire et al. 1996, Martini et al. 1997, Fuente et 
al. 1999, this work) . High densities (n = 10* cm~^) are also de- 
rived from the observations of molecular lines towards the NW 
PDR (Fuente et al. 1993, 1996). As commented in Section 3.2, 
this PDR is formed by high density filaments immersed in a 
more diffuse medium. The observed 01(63 jum)/CII(157.7/xm) 
ratio is very likely determined by the value of this ratio in the 
diffuse component, u ^ 10* cm^"*, which fills most of the LWS 
aperture. In fact, the better angular resolution [CII] 157.7 /im 
and [01] 63 /im data reported by Chokshi et al. (1988) and 
Gerin et al. (1998) show that the emission of both lines is ex- 
tended over the whole HI region. The high density filaments (n 
> 10"' cm~^) are better traced by the observations of the H2 
rotational and ro-vibrational lines towards the NW and SW 
PDRs, as well as by the observation of the rotational lines of 
molecules with high dipole moment. 

i.l. Sill 

The emission of the [Sill] 34.8 /xm line is considered to arise 
mainly in the neutral part of the cloud, i.e., in the PDR. Silicon 
has a first ionization potential of 8.15 eV and remains singly 
ionized until an optical depth of At,~ 5, a value which is fairly 
independent of n and Go. With a critical density of n = 3.4 
10^ cm~^, the intensity of the [Sill] 34.8 //m line is strongly 
dependent on the uncertain Si abundance in gas phase. Silicon 
is heavily depleted in the ISM with logarithmic depletion val- 
ues {Ssi = log(X|"7Xsi)) ranging from -1.30 to -0.18. Most 
PDR models adopt an Si abundance of 3.6 10^^' {5si= -1.0) 
which is the mean value derived by van Steenberg & ShuU 
(1988) but the actual value in the PDR could differ by an or- 
der of magnitude, and consequently the intensity of the [Sill] 
34.8 /im line. The [Sill] 34.8 /xm line is usually compared with 
the [01] 63 /im line because both lines have similar excitation 
conditions and their total column densities are fairly indepen- 
dent of n and Go. Standard PDR models {5si= -1.0) predict 
SiII(34.8/im)/OI(63/im) intensity ratios ranging from ~ 0.06 
- 0.08 for Go = 10^ - 10* and densities ~ 10* - 10® cm"^ (e.g.. 
Burton et al. 1992). 

As it is shown in Fig. 2 and 3, the spatial distribution of 
the Sill emission in NGC 7023 is different from that of all 
the other PDR tracers. All the PDR tracers peak towards the 
vibrationally excited H2 filaments located in the walls of the 
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Table 2. Intensitites of the 01 and CII lines 





NW PDR 


star 


SW PDR 


( X 10 * erg s ^ cm ^ sr ^) 








OI(63Mm) 


7.5±0.1 


3.9±0.8 


3.8±1.3 


01(145.6 fim) 


0.9±0.1 


0.4±0.1 


0.3±0.1 


CII(157.7 ^J,m) 


2.8±0.1 


1.8±0.1 


1.7±0.1 


01(63 Mm)/OI(145.6 fim) 


8.3 


9.0 


12.7 


01(63 Atm)/Cn(157.7 ^im) 


2.7 


2.2 


2.3 



Table 3. SWS observations 



Line 


A( fim) 


NW PDR 


NWl 


Star 


SWl 


SW PDR 


(x 10~^ erg s"^ cm"^ sr"^) 












Recombination lines 


Bra 


4.051 


<2.3 


<1.0 


141.7±1.5 


<0.7 


<1.0 


H2 rotational lines 


H2 S(5) 


6.909 


26.4±3.0 


<4.8 


<7.1 


<3.3 


9.2±3.3 


H2 S(4) 


8.025 


15.0±0.2 


<3.9 


<5.0 


<1.4 


5.6±1.1 


H2 S(3) 


9.665 


40.8±1.3 


<4.4 


<3.3 


<1.6 


21.2±1.5 


H2 S(2) 


12.279 


24.0±3.3 


<10 


<7.4 


<8.0 


16.5±2.7 


H2 S(l) 


17.035 


21.4±0.7 


1.7±0.4 


1.8±0.2 


1.6±0.4 


10.0±0.7 


H2 S(0) 


28.219 


4.8±4.8 


<4.4 


<3.5 


<3.3 


<4.5 


Atomic Fine Structure Lines 


SIV 


10.510 


<2.3 


<2.9 


<3.1 


< 1.8 


<1.7 


Fell 


25.988 


<2.3 


<2.0 


<1.4 


<2.1 


<2.2 


SIII 


33.480 


<3.3 


<6.6 


<8.3 


<1.7 


<3.3 


Sill 


34.814 


3.4±1.3 


4.7±0.5 


8.6±0.1 


1.8±0.3 


<3.4 



cavity in which the star is immersed (the NW and SW PDRs) 
while the Sill emission fills this cavity and peaks towards the 
star. The intense H2 filaments observed towards the NW and 
SW PDRs show that the extinction from the star to these 
positions is lower than 2 mag. Then, the Sill emission is arising 
in the gas located at an extinction of < 2 mag from the star. 

We have compared the Sill data with the OI observations 
reported by Chokshi et al. (1988) which have a FWHM aper- 
ture of 33". Chokshi et al. (1988) measured an intensity of 
1.2 ± 0.1 10"^ erg s"^ cm^ sr"^ towards the NW PDR. The 
SiII(34.8 jj,m) / 01(63 ^m) intensity ratio is 0.03 at this position. 
This value is in agreement with PDR models if the assumed 
Si elemental abundance is ~ 1.5 10~^ {5si= -1.3), i.e., silicon 
is heavily depleted with only the 5% in gas phase . Chokshi 
et al. (1988) observed a position shifted 20" W from the star 
and measured an intensity of 6.2±0.1 10~* erg cm~^ sr~^ . 
Taking into account the dust and gas distribution around the 
star (see Fig. 1), this value can be considered an upper limit 
to the actual 01 emission and we can conclude that the Sill 
(34.8/im)/OI (63 ^m) line ratio is > 0.14 towards the star, i.e. 
a factor of 4 larger than towards the NW PDR. If the Sill 
emission arise mainly in neutral gas, the Si abundance in gas 
phase must be larger by more than a factor of 4 towards the 
star than towards the NW PDR, i.e. at least 20% of the silicon 
is in gas phase towards the star. 

Even more suggestive is to compare the Sill with HI 
emission. The Sill/HI ratio presents a symmetric distribution 
around the star, in which the Sill/HI ratio increases towards 
the star and is a factor of 7 larger towards the star than to- 
wards the NW and SW PDRs. If the Sill arise in the mainly 



neutral layer traced by HI, more than 30% of the Si is in gas 
phase towards the star. Other possibility is that an impor- 
tant fraction of the Sill emission arise in ionized gas which is 
not traced by the HI emission. The contribution of the ionized 
wind to the intensity of the Sill line is expected to be negli- 
gible. We have isarried out some model calculations using the 
code CLOUDYD(Ferland, 1996) assuming the Kurucz stellar 
atmosphere for a star with T^ff — 17000 K, log g=4 cm~^ 
(Kurucz 1979) and the wind parameters derived by Nisini et 
al. (1995). These calculations show that even assuming the so- 
lar abundance for silicon, the Sill intensity in the ionized wind 
is expected to be a factor of 10~^ lower than the intensity of 
the Bra line. Alternatively, the Sill emission could arise in a 
small HII region created by the star. The Sill intensity in an 
ionized sphere around the star with constant density, n — 10^ 
cm~^, could be as high as 3.2 10^^ erg s^^ cm^^ sr^^ if all the 
silicon is in gas phase, i.e., « 40% of the Sill emission towards 
the star. In this case only the 60% of the emission observed 
towards the star would arise in the neutral layer traced by HI. 
But even in this case the Si abundance in the neutral gas must 
increase by a factor of 4 to explain the Sill/HI ratio. 

Then we can conclude that silicon is heavily depleted for 
A„> 2 mag {5si= -1.3) and its abundance increases progres- 
sively towards the star where at least 20% - 30% of the silicon 
must be in gas phase. 

At this point, it is interesting to consider the case of the 
[Fell] 26.0 /im line. We have not detected the emission of the 
[Fell] 26.0 /im line towards any position. This lack of emission 

^We have used the version number 90.04 which uses the MICE 
interface developed by Henrik Spoon at MPE 
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can be explained by excitation effects and a different deple- 
tion of Fc and Si on grains. The critical density of the [Fell] 
26.0 /im line is an order of magnitude larger than that of the 
[Sill] 34.8 //m line (HoUenbach & Mc Kee 1979). On the other 
hand, Si is more easily released to the gas phase than Fe. Fe is 
expected to be significantly depleted in the HII region and the 
FDR although the Si has boon released to the gas phase (Sofia 
et al. 1994). In standard FDR models, the depletion of Fe is 
assumed to be a factor of 10 larger than that of Si {SFe= -2.0) 
and the intensity of the FeII(26.0 /im) lines is then an order of 
magnitude lower than the intensity of the Sill (34.8 /xm) line 
(Burton et al. 1992). 

It is well known that the abundance of gas phase sili- 
con in the interstellar medium increases by grain-grain colli- 
sions and/or sputtering in shocks (Martin-Pintado et al. 1992, 
Caselli et al. 1997, Bachiller & Perez-Gutierrez 1997). Pho- 
todesorption of Si in grain mantles by UV radiation has re- 
cently been proposed as a mechanism to explain the large SiO 
abundance in some PDRs (Walmsley et al. 1999). Our data put 
some constraints on the efficiency of photodesorption to release 
Si to the gas phase. In NGC 7023, Si is heavily depleted (i5s:i= 
-1.3) in gas at an extinction of 2 mag, i.e., silicon is mainly in 
solid form in the molecular gas where the emission of SiO is 
expected to arise. The large bipolar cavity associated with HD 
200775 proves the existence of an energetic bipolar outflow in 
a previous stage of the stellar evolution (Fuente et al. 1998). 
The large gas phase Si abundance in the positions closest to 
the star could be a relique of this stage. If it were the case, the 
efficiency of photodesorption to return the Si to the gas phase 
would be even lower. 

4-2. A gradient in the ortho-to-para-H2 ratio 

The S(0),S(1),S(2), S(3), S(4) and S(5) v=0-0 Ha rotational 
lines have been detected towards the NW PDR and the 
S(1),S(2), S(3), S(4) and S(5) lines towards the SW PDR. In a 
previous paper (Fuente et al. 1999), we have studied in detail 
the H2 rotational lines towards the NW PDR. In this paper we 
will concentrate in the observations towards the SW PDR. 

In Fig. 5 we show the rotational diagram of the H2 lines 
towards the SW PDR. Note that the errors in Fig. 5 are en- 
tirely dominated by the calibration uncertainties (15% at 6.909 
urn, 25% at 8.025 /im, 25% at 9.665 fim, 25% at 12.279 %, 
20% at 17.035 /im and 30% at 28.219 /xm). The data have 
been corrected for dust attenuation using the value for the 
dust opacity derived from the LWSOl spectra (A„ = 0.3 mag) 
and the extinction curve of Draine & Lee (1984). The extinc- 
tion for the S(0), S(l), S(2), S(3), S(4) and S(5) lines amounts 
to 0.006,0.013, 0.013, 0.027,0.013 and 0.006. Like in the case 
of the NW PDR (see Fuente et al. 1999), the rotational di- 
agram of the SW PDR shows that the ortho-H2 levels have 
systematically lower Nu/fju values (where Nu and <;„ are the 
column densities and degeneracies of the upper levels of the 
transitions) than the adjacent J-1 and J-|-l para-Ha levels pro- 
ducing a "zig-zag" distribution. In fact, the ortho- levels seem 
to define a curve which is offset from that of the para-levels by 
more than a factor of 2 (see Fig. 5). This offset cannot be ex- 
plained by calibration uncertainties and extinction effects. The 
calibration uncertainties are at most 30%. Since the values of 
the extinction at 17.03 ^tm, 12.3 /im and 8.025 /xm are very 
similar (Draine & Lee 1984), a correction for dust attenuation 
cannot push the S(l), S(2) and S(4) points to the same curve. 



This " zig-zag" distribution cannot be due to the different aper- 
tures of the instrument for the different lines. In the case of a 
point-like source, the S(l) and S(2) lines should be corrected 
by a factor 1.35 relative to the S(3) line and the offset between 
ortho- and para- curve would increase. 

The "zig-zag" distribution observed in the rotational dia- 
gram can only be fitted by assuming a non-equilibrium OTP 
ratio, i.e., different from the OTP equilibrium value at the gas 
kinetic temperature. In this case, only the rotation temper- 
ature between the levels of the same symmetry constitute an 
estimate of the gas kinetic temperature. Using the ortho-levels, 
we have derived a rotation temperature of ~ 450 K between 
the S(l) and S(3) lines and of ~ 630 K between the S(3) and 
S(5) lines. For the para- levels we have derived a rotation tem- 
perature of ~ 430 K using the S(2) and S(4) lines. Then we can 
conclude that the H2 rotational lines are arising in gas with ki- 
netic temperatures ranging from ~ 400 to 700 K. These kinetic 
temperatures are similar to those obtained in the NW PDR, 
but the intensities of the lines are about a factor of 2 lower. A 
change in the incident UV field and/or the density would imply 
a variation of both the line ratios (rotation temperatures) and 
the line intensities. Since the line ratios are similar to those of 
the NW PDR and the line intensities are significantly lower, 
we conclude that the SW PDR has the same excitation condi- 
tions as the NW PDR (Go = lO", n = lO'' cm'^) but a beam 
filling factor of ~ 0.5. All PDR models which include both, 
chemical and thermal balance in the gas, assume a fixed value 
of 3 for the OTP ratio, which is the equilibrium value for gas 
kinetic temperatures > 100 K. In order to determine the value 
of the OTP ratio in the SW PDR, we have corrected the in- 
tensities predicted by the model of Burton et al. (1992) for 
the effect of different values of the OTP ratio assuming that 
the total amount of H2 molecules at each temperature and the 
line ratios between the levels of the same symmetry remain 
unchanged. This correction is only valid if the excitation of 
the H2 lines is mainly coUisional which is the expected case 
for the low rotational transitions. In Fig. 6, we compare the 
observational data for the S(0), S(l), S(2) and S(3) transitions 
with the model for Go = lO", n = lO" cm"-^ and a fixed OTP 
ratio of 3, and the same model corrected for OTP ratios of 1.5 
and 2. We obtain that the rotation diagram of the SW PDR 
is well fitted with an OTP ratio of 1.5 (see Fig. 6). Martini et 
al. (1997) derived an OTP ratio of 2.3±0.5 towards the SW 
PDR based on the H2 vibrational lines. Because of the opti- 
cal depth effects in the excitation of the ortho- and para- H2 
vibrational lines, this value is just a lower limit to the actual 
value of the OTP ratio (Draine & Bertoldi 1996, Sternberg & 
Neufeld 1999). Then we can conclude that the OTP ratio is 
larger than 2.3±0.5 in the less shielded layers of the PDR (A„ 
< 0.7 mag) where the vibrational lines arises. The gas kinetic 
temperature in these layers ranges from several hundreds to 
2000 K. Comparing the OTP ratio derived from the rotational 
lines with that derived from the vibrational lines, we conclude 
that the OTP ratio increases across the PDR. It takes values 
^ 1.5 at temperatures of ~ 400 - 700 K and values close to 3 
at temperatures ~ 700 K. 

The same behavior, a non-equilibrium OTP ratio for the 
gas with kinetic temperatures ~ 400 K and a gradient of the 
OTP ratio across the PDR was found in the NW PDR by 
Fuente et al. (1999). They discussed that the most likely expla- 
nation for this behavior is the existence of an advancing pho- 
todissociation front. An important problem remains in this ex- 



6 



planation, the velocity of the photodissociation front required 
to have a non-cquiUbrium OTP ratio is too large (~ 10^n~^ 
kms^^). However, since the expected behavior of the OTP ratio 
in an advancing photodissociation front is in agreement with 
that observed in our data (the cooler gas is further from the 
OTP equilibrium value than the warm gas) and there are many 
uncertainties in this estimate (it is not known the OTP ratio 
at which the H2 is formed within the PDR, the gas is clumpy 
and small dense clumps can be evaporating into the PDR,...), 
we considered that this was the most likely explanation. Later, 
Lemaire et al. (1999) have obtained a very high spectral resolu- 
tion spectrum of the H2 S(l) v=l-^0 line at 2.121 /im towards 
the NW PDR. The H2 line is seen shifted in velocity relative to 
that of the HI and the CO rotational lines. They argue that this 
velocity difference is a clear proof of the presence of dynamical 
effects in this PDR and estimated a (projected) velocity of 1 
kms~^ for the photodissociation front. So far, very few mod- 
els have been developed for non-equilibrium PDRs (Bertoldi & 
Draine 1996, Storzer & HoUenbach 1998). These models con- 
sider regions illuminated by O stars in which a ionization front 
is combined with a photodissociation front and assume a fixed 
value for the OTP ratio of 3. They do not give any informal 
tion about the effect of an advancing photodissociation on the 
ortho-to-para- H2 ratio. But they predict how the existence of 
a photodissociation front affects the thermal structure of the 
region. Following the model by Storzer & HoUenbach (1998), 
because of the advection of molecular gas through the PDR, 
the intensities of the H2 rotational lines and vibrational lines 
can be enhanced by a factor of 3. Assuming a velocity of 1 
kms~^ for the photodissociation front, the intensities of the 
H2 rotational lines observed in the NW PDR can be well fitted 
with a density of n- lO'^ cm"^ and a UV field of Go = lO"*. A 
non-equilibrium OTP ratio has to be assumed to explain the 
"zig-zag" distribution. 

4-.3. Absence of molecular emission 

Towards NW PDR we have tentatively detected a weak line 
centered at 153.927 /xm with a S/N ratio of ~ 3. A weak line 
centered at 153.637 fjxn appeaxs also in the LWS spectrum 
towards the star position with a S/N ratio of 5. Taking into 
account that the spectral resolution of the LWS at these fre- 
quencies is 0.6 fim, both lines can be identified with the CO 
J = 17— >16 line (153.267 /im). However, since other CO lines 
have not been detected towards these positions, we cannot dis- 
card the possibility of a misidentification or a spurious detec- 
tion. If confirmed, this line would implied the existence of a 
gas with Tfe> 200 K and n> 10'' cm^'^ in the viciruty of the 
star. However the density of the bulk of the molecular gas in 
the PDR must be at legist an order of magnitude lower. Fe- 
derman et al. (1997) estimated based on the ultraviolet ab- 
sorption lines of CH, CH"*", C2 and CN that the density of 
the PDR in front of the star is ~ 200 cm~^. Furthermore, the 
01(63 /xm)/CII(157.7/xm) ratio is consistent with densities of 
10"* cm-3. 

We have not detected any OH, CH and CH2 line towards 
the observed positions. Emission of the OH lines have been 
detected in several sources (NGC 7027: Liu et al. 1996; R CrA 
and LkHa 234: Giannini et al. 1999). ISO-LWS observations 
towards a sample of 11 Herbig Ae/Be stars shows that the 
OH lines are only detected towards the stars associated with 
with the largest mean densities (n> 10® cm~^) as derived from 



the 01(63 ^im)/CII(157^m) ratio. The non-detection of OH 
in this nebula is consistent with the idea that the bulk of the 
molecular gas is formed by the low density component traced 
by the OI(63 Atm)/CII( 157.7 ^tm) ratio. 

A similar situation is found for CH and CH2. Large CH 
column densities have been detected in absorption. The data 
by Federman et al. (1997) shows that the CH column density 
in front of HD 200775 is 3.2±0.2 10^^ cm"^. Even ff we assume 
excitation temperatures as low as Tex = 20 K, we should have 
detected intensities of the order of 10~® erg s~^ cm~^ for the 
180.93 and 149.09 /xm lines. The non-detection of these lines 
can only be explained invoking the excitation conditions. The 
beam filling factor of the dense gas (n> 10'' cm^"^) is very low. 
Another possibility is that the CH abundance decreases for 
high densities. This would also explain the lack of CH emission 
in other sources. In fact, far to our knowledge, emission of CH 
has not been detected towards any source (just a tentative 
detection towards NGC 7027). 

5. Summary: A model for the nebula 

In this paper we present the ISO-SWS and LWS observations 
towards a strip crossing the star and the most intense PDRs in 
the reflection nebula NGC 7023. These observations altogether 
with previous maps in molecular lines and the [HI] 21cm line 
have allowed to figure out a model of the nebula. 

We can consider the nebula as composed by three different 
regions. The first one is formed by the star, the small HII region 
around it and the low density gas filling the cavity. Federman 
et al. (1997) based on studies of absorption lines derived that 
the density of the gas in front of the star is of a few 10^ cm~^. 
A slightly larger value, n — 10^ cm~^, was found by Gerin et al. 
(1998) based on observations of the CO rotational lines towards 
the star. The Bra and [Sill] 34.8 /im lines peak in this region 
but the [HI] 21cm and H2 rotational lines have a minimum in 
emission. The Bra emission is well explained as arising in the 
stellar ionized wind modeled by Nisini et al. (1995). However, 
the contribution of the stellar wind to the intensity of the [Sill] 
34.8 fim line is expected to be negligible. The [Sill] 34.8 /im 
line is better explained as arising in the diffuse gas filling the 
cavity in which at least 20% - 30% of the silicon is in gas phase. 

The second region is formed by the walls of the cavity in 
which the star is immersed. Clumpy PDRs have been formed in 
these walls with high density filaments (n ~ 10'' - 10^ cm^"*) 
immersed in a more diffuse medium( n ~ 10* cm~'^). These 
high density filaments are located almost symmetrically ~ 50" 
NW and ~ 70" SW (NW PDR and SW PDR) from the star 
and constitute the peaks of the [CII] 157.7 fim, [01] 63.2 and 
145.6 fiin, [HI] 21cm and H2 rotational lines emission. The NW 
and SW PDRs have very similar excitation conditions although 
the SW PDR is a factor of 2 weaker than the NW PDR. This 
difference in intensity is very likely due to a different beam 
filling factor. In both, the NW and SW PDRs, the intensities of 
the H2 rotational lines can only be fitted assuming an ortho-to- 
para-H2 ratio lower than 3. Since the rotation temperatures of 
the H2 lines ranges from 400 to 700 K, this implies the existence 
of a non-equilibrium OTP ratio in the region. The comparison 
between the vibrational and pure rotational H2 lines suggests 
that the OTP ratio increases across the PDR from values as 
low as ~ 1.5 for gas with temperatures ~ 400 - 700 K to values 
close to 3 in the gas located at a visual extinction A„ < 0.7 
mag {Tk ~ 700 K). This non-equilibrium OTP ratio has been 
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interpreted in terms of an advancing photodissociation front. 
The intensity of the [Sill] 34.8 fim line observed towards the 
NW PDR is consistent with the predictions of PDR models 
with a standard silicon depletion of 5 Si = -1.3, i.e., only 5% 
of the silicon is in gas phase in these PDRs. 

The third region is the molecular gas. We have not detected 
the OH, CH and CH2 molecular lines towards this source. This 
is consistent with the weakness of these lines in other sources 
and reveals that the beam filling factor of the dense gas (n > 
10^ cm-3) in the LWS aperture is low. The CO J=17^16 fine 
has tentatively been detected towards the star. 
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Fig. 1. Scheme of the observed region. The HI total column density map (contours) is superposed on the integrated intensity map of the 
^■^CO J=1^0 map (grey scale) (Fuente et al. 1998). The filaments observed in H2 fluorescent emission are also shown (Lemaire et al. 1996). 
The observed positions axe indicated by filled squares and the star. 
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Fig. 2. Spectra of the H2 S(0),S(1),S(2),S(3),S(4) and S(5) rotational lines , the Sill (34.8 fim) line and the Bra HI recombination Une 
towards the five positions observed across the PDR associated with NGC 7023 (see Fig. 1). 
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Fig. 3. Normalized intensity of H2 S(l), HI 21cm and SiII(34.8 /im) lines as a function of the distance from the star. Note that both, the 
H2 and HI emission peaks towards the H2 filaments, while the SiII(34.8 /im) emission peaJfs towards the star. 
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Fig. 4. LWS spectra towards the NW PDR and the star position after subtracting the continuum emission. 
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Fig. 5. Rotational diagram of the pure H2 rotational lines observed towards the SW PDR. The points corresponding to para- levels are 
joined by a dashed line, and the points corresponding to ortho-levels with a thick line. Note that there is a shift between the ortho- and 
para-H2 lines that is larger that the observational errors. 
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Fig. 6. Comparison of the observational data with the results of PDR models with Go = 10*, n = 10^ cm^"^, and OTP ratios of 3 (open 
triangles), 1.5 (open circles) and 1 (open squares). The data are shown with black squares and errorbars. The beam filling factor has been 
taken 0.5 for all the models. The model with OTP = 1.5 fits the observational data. 
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